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A difficult problem in quantum network communications is how to efficiently transmit quantum 
information over large-scale networks with common channels. We propose a solution by developing a 
quantum encoding approach. Different quantum states are encoded into a coherent superposition state 
using quantum linear optics. The transmission congestion in the common channel may be avoided by 
transmitting the superposition state. For further decoding and continued transmission, special phase 
transformations are applied to incoming quantum states using phase shifters such that decoders can 
distinguish outgoing quantum states. These phase shifters may be precisely controlled using classical chaos 
synchronization via additional classical channels. Based on this design and the reduction of multiple-source 
network under the assumption of restricted maximum-flow, the optimal scheme is proposed for specially 
quantized multiple-source network. In comparison with previous schemes, our scheme can greatly increase 
the transmission efficiency. 

Maximizing information transmission is very important concern when dealing with limited- resource 
scenarios in large-scale networks. Network coding', which allows multiple messages to be encoded 
before transmission in common channels, may provide a solution for single-source networks. One 
typical example is illustrated in Figure 1(a). With linear encoding, the source node can simultaneously transmit 
two messages to all receivers from two edge-disjoint paths', even if two receivers share one common channel CD. 
The key is that two incoming messages {x, y} are encoded into a new message x + y,zt the node C. Each receiver 
can recover {x, y} using {x, x + y} [y, x + y). This task cannot be successfully performed using the trivial 
transmission scheme [store-forward routing for each node^] because in this scheme, only one message, x ory, may 
be transmitted in the common channel CD each time. Generally, transmission conflicts in common channels 
[network congestion] become serious problems [congestion collapse] in large-scale networks [Internet]^. 
Fortunately, network coding can achieve the optimal Shannon capacity of single-source networks''^. Network 
coding is considered to be an important technology for next-generation communication to achieve network 
multicasts''"'' [the sender simultaneously sends multiple messages to all receivers on a single-source network] or k- 
pair transmissions^ ", i.e., multiple unitcasts [each sender transmits messages to its corresponding receiver simul- 
taneously over multiple-source networks]. 

When classical networks are quantized with quantum nodes and quantum channels, one may wish to optimize 
the transmission efficiency over large-scale quantum networks, as illustrated in Figure 1(b). Because of the 
quantum no-cloning theorem, the source node O may be replaced with two source nodes Si and S2. Thus, the 
multicast over the Butterfly network is reduced to the 2-pair problem on the reduced quantum network with 
common channels'. However, perfect quantum A:-pair transmissions have been proved to be impossible in the 
absence of any additional resources'""'^. They require all incoming states |(^,) be encoded into a coherent 
superposition at the node C, and decoded at the node D. This is a rather difficult task. The situation 

may be changed if classical communications are freely allowed'''''''. This assumption is reasonable because 
classical communications are much cheaper and more readily available than quantum communications. These 
results are dependent on the linear encodings''' or the nonlinear encodings'^ applied in the enlarged encoding 
space, and give rise to three natural questions. The first question (Ql) is whether the enlarged encoding space is 
necessary for large-scale quantum network communications. The second question (Q2) is what amount of 
classical communication is sufficient. Classical two-way unlimited channels have been assumed to exist between 
any two nodes''' or only between two quantum nodes connected by quantum channels'^. The third question (Q3) 
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Figure 1 | Schematic network transmission over the Buttery network, (a) 

The classical network multicast with network coding. Each channel has 
unit transmission capacity. Each node _R, has two edge-disjoint paths 
originating from the source node 0,i= 1,2. CD is a common channel for 
two receivers. The input messages {x, y] are encoded into x + y at the node 
C, and forwarded to the node D. x + y is copied [unit fidelity], and each 
message is forwarded to one node R;. Each receiver can recover {x, y] from 
{x,x + y] or {y, x + y}. (b) The quantized network without the node O. All 
nodes and channels may be quantized with quantum participants and 
quantum channels, respectively. The node O is canceled because of the 
quantum no-cloning theorem at the node S,. Thus, the quantum task is for 
each Si to multicast an unknown state to all receivers simultaneously, (c) 
The reduced quantized network. The transmissions over the channels S1R2 
and S2-R1 are trivial and reduced. The remaining task is to design the 
transmission over the common channel CD. 

is whether the linear encoding is sufficient for large-scale quantum 
network communications. These problems are related to the optim- 
ization of the transmission capacity of multiple-source quantum 
networks. 

In this paper, based on the achievements of quantum information 
theory"'"^" and quantum networking theory^' we investigate these 
problems using similar ideas in classical network coding. Unlike the 
entanglement quantization of a classical channel'"", it may be quan- 
tized with a continuous physical channeP*"', and the transmission 
information may be quantized with quantized electromagnetic fields 
of identical frequencies. To distinguish the decoded quantum states 
for different nodes, special phase-shift operations may be designed to 
index different incoming quantum states, using phase shifters [coup- 
ling the optical fields to driven Duffing oscillators] or gauge trans- 
formations'". This encoding can spread the spectral content of the 
quantum information across the entire spectrum in order to encode 
the information, and can distinguish different senders with their own 
phase factors. Unfortunately, the added phase information is not 
easily decoded, and doing so requires the ability to precisely control 
the remote chaotic systems during the communication. To solve this 
problem, free classical communications are assumed between two 
quantum nodes with common channels, and are used to synchronize 



the remote phase shifters. This is classical chaos synchronization and 
may be achieved using the nonlinear coupling between the optical 
fields and Duffing oscillators" or semiconductor lasers'" The 
key is the nonlinear Kerr interaction, which can be used to couple the 
classical chaotic light with the information-bearing quantum lighf^. 
Recently, electrooptic modulators (EOMs) have also been used for 
chaos synchronization""" Moreover, all encoded quantum 
information may be combined into a coherent superposition state, 
and decoded using paired multiport beam splitters'*'^"'''. Based on our 
transmission schemes in common channels, using the network 
reduction shown in the supplementary information (SI), we can 
identify the optimal transmission scheme for quantum multiple- 
source networks assuming restricted maximum-flow. Our scheme 
is beyond both the quantum fc-pair transmissions" based on the 
classical solvability and classical network transmissions^^"^'' via net- 
work coding. These results maybe beneficial for large-scale quantum 
network communications. 

Results 

Consider an acyclic directed quantized network Gq= (Vq,!?^), as 
shown in Figure 2(a)). and £q are the node set and the edge 
set, respectively. Each node in Vq is quantized with a quantum 
participant that can perform all quantum operations and classical 
operations. The transmission information is quantized with electro- 
magnetic fields Uj and bj of the same frequency. Each channel in £q is 
quantized with a continuous-time physical channel and has a unit 
transmission rate [one quantum state]. Each pair (S„ Rj) has Ijj a 1 
edge-disjoint paths. Our task is to allow efficient transmissions in 
large-scale quantum networks: All source-sink nodes pairs commun- 
icate simultaneously, subject to restricted maximum-flow. Although 
quantum multiple-source networks have no uniform network topo- 
logy, based on the network reduction shown in the SI, the optimal 
quantum multiple-source transmission may be reduced to the trans- 
mission in the primitive network, as shown in Figure 2(b). Thus, 
special encoding and decoding operations should be designed to be 
suitable to the quantum task. 

Restricted maximum-flow. In the optimization theory, the maximum- 
flow problem [unit capacity for each channel] is equivalent to 
identifying the maximal number of edge-disjoint paths between the 
source and the sink, under the assumption unit capacity of per edge. 
Our interest is in the case of restricted maximum-flow, i.e., no 
common channels for different source-sink node pairs are outgoing 
edges of source nodes or incoming edges of sink nodes. This 
assumption is reasonable because of the quantum non-cloning 
theorem. 

Motivated by network coding' and quantum network theory^'"^^, 
a schematic illustration of quantum transmission over a common 
channel is presented in Figure 3. The information-bearing field Uj 
originating from the node A, is first shifted at the node C using a 
chaotic phase shifter (CPS^) with the Hamiltonian 5j{t)ajaj and the 
time-dependent classical chaotic signal 5j{t),j= 1, • ■ ■ ,n. This phase 
shift corresponds to the gauge transformation in the nearest nodes^". 
All new quantum information is encoded using a multiport beam 
splitter (MBSi), and transmitted over the common channel CD. The 
combined quantum information is amplified using a phase-insens- 
itive linear amplifier (LA) at the node D to compensate for the 
information losses induced by MBSi, and then decomposed into n 
different components by MBS2. The amplifier gain of the LA is n + L 
All decomposed information may be decoded using CPS'; [the 
inverse of CPS^] with the corresponding Hamiltonian — 5'-{t)bjbj, 
j= 1, ■ ■ ■ ,n. Each decoded information-bearing field bj is sent to the 
subsequent node Bp j = I, ■ ■ ■ ,B„. 

Note that each pair of CPSj and CPS'^ induces phase shifts with the 
phase exp{iOj{t)) and the inverse phase exp{ — iOj{t)), respectively. 
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Figure 2 | Schematic acyclic directed quantum multiple-source network. 

(a) {Si, • • • Sk}, {Ri, ■ ■ ■ ,-R(} t= V are source and target nodes, respectively. 
Each pair SiRj has lij edge-disjoint paths. CD is a common channel for 
different pairs. No common channel is the outgoing edge of a source node 
or the incoming edge of a sink node, (b) The primitive subnetwork of 
quantum multiple- source networks, aj and bj are information-bearing 
fields of quantum information [the original fields generated by the former 
nodes] . The node C has n incoming edges, and the node D has n outgoing 
edges. The transmission in the channel CD is the main concern for 
quantum multiple- source transmission tasks. All incoming states 
l*^!)) ■ ■ ■ Mil) should be encoded into a superposition state at the node C, 
and decoded at the node D. 



where Oj{t)= I 5j{T)dz. To achieve faithful transmission, these 
J 0 

transformations must be precisely controlled during the process of 
quantum communication, i.e., it must be ensured that the two cha- 
otic systems have the same parameters, initial values, and evolutions 
such that 5j{t) = 5'j{t) for each j= 1, ■ ■ ■ ,«. However, this precise 
control is impractical for remote participants because the chaotic 
system is unstable in system parameters and initial values. 
Therefore, additional classical channels are assumed to exist for 
common channel CD, and used to synchronize each pair of CPSj 
and CPSj,; = l,--- as shown in Figure 3(b). These classical chan- 
nels are cheap and readily available compared with quantum 
channels. 

Modeling quantum transmission over a common channel. Con- 
sider the primitive quantum network shown in Figure 3(b). Each pair 
of CPSj and CPSj has been synchronized prior to the transmission of 
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Figure 3 | Schematic quantum transmission over a common channel, (a) 

Quantum transmission without additional classical channels. fl|^ is the 
creation operator of the auxiliary vacuum field entering the LA, and 
denotes the annihilation operator entering the second MBS (MBS2), 
j= 1, • • • ,M. The encoding at the node Cis performed using CPSyandMBSi. 
The decoding at the node D is performed using the LA, MBS2 and CPSj 
[the inverse of CPSy]. (b) Quantum transmission with additional classical 
channels for the common channel. The encoding is identical to that shown 
in Figure 3(a)), whereas the decoding is performed using the chaotic 
synchronization of CPSj and CPSj, j= 1, • • • 

quantum information, j=l, ■ ■ ■ ,n. The information-bearing fields 
ai,---,a„ with the same frequency oj^ are modulated using n 
different pseudo-noise signals and pass through CPS^, MBSi, the 
LA, MBS2, and CPS'j sequentially. The global quantum transmis- 
sion can be described using the following linear relation: 



for all j=l, - ■ ■ ,n. The matrices {ciij)„^„ and (pij)nxn denote the 
transformations of MBSi and MBS2, respectively, and satisfy 
/lau/iii — ■ ■■ =AaLi„P„^ = 1. fls3 denotes the annihilation operator 
of the auxiliary vacuum field entering MBS2, s = 2, ■ ■ ■ For the 
pseudo-noise chaotic phase shift Oj{i) from the CPSs, one needs to 
take an average over the broadband random signal, i.e.. 



exp( 



my 



with Mj = exp 



- n I Sa (co) / m doj 
' / / 

and the power-spectrum density Ss^{oj) of signal Sj{t). Here, coij 
and ojjy are the lower and upper frequency-band bounds of Sj{t), 
respectively. Thus, the equation (1) is further reduced to 



hi = a: 



1 



J2 sjMjMkak + J2 V^jPp'^BS.-^ 



1 



All Mj are extremely small with respect to the chaotic signal with the 
broadband frequency spectrum, and thus can be ignored in equation 
(2), therefore 



(3) 
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i.e., faithful transmission of quantum information is achieved from 
the node Aj to the node A^', j= 1, ■ ■ ■ 



Quantum state transmission over a common channel. Consider 
pure qudit state transmission over the proposed model, as shown in 
Figure 4. The transmitted states are dark states of general A-type d + 

1-level atoms j^^^ = ^^I'^J ^''^ ^ -'^l' where the j-th 
atom is located in the cavity CAj [see Figure 4(a)], j= 1, ■ ■ ■ These 
states are transferred to the cavities via Raman transitions, 
transmitted over the quantum network, and stored in the new 
cavities. Assume that In coupled atom-cavity systems have the 
same parameters. By adiabatically eliminating the highest energy 
level I the atom j will always lie in the lowest d energy levels 
|0)^, ■ ■ ■ — 1)^. The neighboring transition >|i-|- 1)^- is driven 
by a near-resonant laser field, and is coupled to the classical control 
field and the quantized cavity field with a coupling strength Qij(t). 
The Hamiltonian of the atom-cavity systems can be expressed as 



d-i 



(4) 



where Cj is the annihilation operator of the j-th cavity mode; gjj{t) = 
gQij{t)/A is the coupling strength tuned by the classical control field 



-level atom, with d neighboring transitions 




Figure 4 | Quantum state transmission over a common cliannel. (a) The 

circular symbols denote atom cavities. Q.ij( t) denotes the amplitudes of 
classical driving fields in each cavity. Clij denotes the transition frequency 
of \ to \ i + The circle in the center denotes a general A-type d + 1- 
level atom, ,„ [vacuum states] is the input field, j= (b) 
Schematic quantum transmission over a common channel. Uj ,„ denotes the 
incoming information-bearing fields of nodes Aj,j=l, ■ ■ ■ ,n. aji and 13^2 are 
the incoming and outgoing information-bearing fields of MBSi, 
respectively; and Uj^ and Uj^ are the incoming and outgoing information- 
bearing fields of MBS2, respectively, j = 1, • • • and flu are one of the 
incoming and outgoing information-bearing fields of the LA, respectively. 



Slij(t), j = 1, ■ ■ ■ ,n; and A is the atom-cavity detuning. Cj is related to 
the traveling field Uj as follows: 



j.out ■ 



-■V^iCj+Ujjn, 

■■VJ^c'j + a'j, 



(5) 
(6) 



where A is the decay rate of the cavity field. 

CPSy and CPS'j are realized by coupling the optical fields to Duffing 
oscillators" as described by the Hamiltonian 



-n, 



(7) 



where pj and qj are the normalized position and momentum of the 
Duffing oscillators, respectively; cuq is the frequency of the fun- 
damental mode; and ^i, )', and oj^ are constants. The interaction 
between the field Uj and Duffing oscillators is given by the 
Hamiltonian Hj — l^jX^ajaj, where is the coupling strength between 
the field Uj and the oscillators. By choosing a suitable interaction, a 

phase factor exp ^ — !^ CjXj{T)dT^ can be generated for the field a^. 

Moreover, the chaotic synchronization between CPS^ and CPSj may 
be achieved by using the harmonic potential coupling 

To show the quantum transmission efficiency, let us calculate the 

fidelity Fj = (^(j>j\^p'^(j>j^, where p'j is the quantum state received by the 

atom j'. From equation (2), the fidelity F, can be approximated as 1 
when M = 0, i.e., when Duffing oscillators enter the hard chaotic 
regimes''^ This result means that qudit states can be faithfully 
transmitted over this primitive network via a common channel. 

Quantum transmission in a multiple-source network. Note that 
according to the network reduction shown in the SI, the number of 
incoming channels for one common channel is equal to the number 
of outgoing channels. Thus, each common channel CD is equivalent 
to m distinct quantum channels QDi, ■ ■ ■ .CmD^ [not common 
channels] aided by additional classical channels, where m denotes 
the number of incoming channels, as shown in Figure 5. By replacing 
all common channels with equivalent quantum channels, an 

equivalent multiple-source network = (v^,£q^ can be 

constructed, which satisfies that all pairs (Sj, Rj) have no common 
channels under the assumption of restricted maximum-flow. Here, 
the source nodes and the sink nodes are unchanged. The resultant 
quantum transmission can be easily achieved via forward routing on 
the equivalent network with the aid of chaotic synchronization on 
the auxiliary classical channels. Thus, we identify and implement the 
optimal quantum transmission under the assumption of restricted 
maximum-flow in a multiple-source network, and partially answer 
the questions Q1-Q3. More specifically, the un-enlarged linear 
encoding [encoding operations such as those represented by 
equation (1)] is sufficient for large-scale quantum network 
communications under the assumption of restricted maximal-flow, 
and unlimited classical one-way channels corresponding to the 
common quantum channel are assumed. Of course, classical 
synchronization should be applied prior to the transmission in the 
common channel and requires some classical communication. 

Discussion 

We have introduced quantum multiple-source networks based on 
classical multiple-source networks and chaotic synchronization, 
where quantum information can be simultaneously transmitted in 
multiple subnetworks derived from source-sink node pairs. The 
proposed quantum transmission attains the optimal transmission 
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Figure 5 | Schematic illustration of quantum transmission in a multiple- 
source quantum network. Si, - ■ ■ ,8^ are source nodes, and Ri, • • • ,R, are 
sink nodes. The channel CD is a common channel for various source-sink 
node pairs. CiDi, • • • ,C„,D,„ are the equivalent m distinct quantum 
channels, supported by additional classical channels, m denotes the 
number of distinct quantum paths. All pairs ( S,-, Rj) have no common 
channels in the equivalent network under the assumption of restricted 
maximum-flow. 

capacity under the assumption of restricted maximum-flow; in this 
respect, it is superior to other proposed approaches' '\ Unlike 
quantum fc-pair tasks with enlarged linear encodings'"' ''' or nonlinear 
encodings'^, our encoding is linear and is not completed in an 
enlarged space. And classical channels are assumed to exist only 
for the common quantum channels and not all quantum channels'^. 
Moreover, the new scheme extends the solvability of the quantum 
fc-pair problem, and is more general than previous schemes" '^ 
which depend on the solvability of the classical A:-pair problem. 
One typical example is presented in the Figure 1(c); this example is 
unsolvable'"^"^'' in terms of the classical 2-pair problem using network 
coding. Other examples are provided in the SI. Our result is also 
different from continuous-time quantum walks'^, these examples 
break the time-reversal symmetry of the unitary dynamics for the 
purpose of enabling directional control, enhancement, and suppres- 
sion of quantum transport. However, it requires controlling the 
chiral system for the gauge transformations, which may be difficult 
in case of remote systems in large-scale quantum networks. We use 
auxiliary classical channels to achieve similar transformations with 
the aid of chaotic synchronization. The proposed scheme can avoid 
to precisely controlling different chaotic phase shifts, and should be 
useful for large-scale quantum networks. Furthermore, the linear 
optical equipments are used to encode different quantum states 
and solve the transmission congestion on common channels. 
GeneraUy, we can get the optimal transmission scheme in mul- 
tiple-source network under the assumption of restricted max- 
imum-flow. Of course, our proposal entails three important 
experimental requirements. The first one is the quantum interference 
of signals from different chaotic sources"'^'^'". The second is the imple- 
mentation of chaotic phase shifters and their synchronization, which 
may be achieved using all-optical systems^'' or optoelectronic''*. The 
third is the implementation of multiport beam splitters''""^'. More 
specifically, only the coefficients of the first column or row are rel- 
evant to our proposal; this feature may reduce the design complexity. 



These issues may not be far out of reach because of recent experi- 
mental and theoretical developments. Our schemes may provide one 
method for long-distance quantum network communications. 

Methods 

We calculate the linear mapping over the quantum network shown in the Figure 4(b). 
The mapping relationships of the CPS'i , ■ ■ ■ ,CPSn may be represented as 

/ "11 \ \ / "' \ 



(8) 



where and a, are the annihilation operators of the auxiliary vacuum fields entering 
and output the CPS, respectively, / — 1, ■ ■ ■ ,n. The mapping relationship of the MBSi 
is defined as 



/ 012 \ 

\a„2 J 



( an 



«-ln\ ( flu \ 

Vfl»i/ 



\ =<nl ■ • • ««» / 

where (3i2, ■ ■ ■ ,fl„2 the annihilation operators of the auxiliary vacuum fields 
entering the MBSi- The mapping relationship of the LA is defined as 



313 = Aai2 + V G — A^ajjj, 
3l3 = V G-X^a\2+WA, 



(9) 

(10) 
(11) 



where a^^ is the creation operator of the auxiliary vacuum field entering the LA. The 
mapping relationship of the MBS2 is defined as 



\an4 / 



(P. 



\fln3 / 



(12) 



where ajj and ^,4 are the annihilation operators of the auxiliary vacuum fields entering 
and outcoming the MBS2 respectively, i— 1, ■ ■ ■ ,n. The mapping relationship of the 
CPS'i, ■ ■ ■ CPS'„ are defmed as 



(13) 



where bj is the annihilation operator of the auxiliary vacuum field outcoming the 
CPS,, j— 1, ■ ■ ■ Then, combining these equations and the conditions ?. = n, 
an — ■ ■ ■ — y.in — Pu— ' ' ' — Pni — l/\/"i the total input-output relationship of the 
quantum network is 



Ai^l 


■(' J 


/ 014 \ 


\b„) 




\fln4 / 



where 9j is independent chaotic noise, j— 1, 



3 -H i/ e-'ai 



(14) 



Averaging the chaotic phase shift. The chaotic signal 3j{t) may be expressed as a 
combination of many high-frequency components, i.e.. 



(15) 



where Aji^, co^p-, are the amplitude, frequency, and phase of each component, 
respectively. Then the phase of the signal is defined as 



Using the Fourier-Bessel series identity^^ j^exp(fx sin y) — nJni^) exp(m)')j 
with the «-th Bessel function of the first kind we can write 



exp(-;e,-(t))=n 



(16) 



Take average over the random phase o}^^, the components related to the high-fre- 
quencies is averaged out because of the energy dissipation. It means that the resultant 
is only the near-resonance components, i.e, the lowest-frequency terms [tt^ — 0] 
dominating the dynamical evolution. Thus, we have 



exp(-,'0,(O)-n/o(^ 



(17) 
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Moreover, since the chaotic signal Si{t) is mainly distributed in the liigh-frequency 
regime, we have Ajk<^cOjk. 



(18) 



Using the approximations Jo{x) ~ 1 — and log(l + x) ~ x for x<^l, it easily 
follows that 



n/o 

k 



' exp 



(19) 



(fo) 



= exp^ 

where Mj — exp ^ — 7c^ 5,^^((y)/ci)^(ioj^ . 

Consequently, from equations (9) and (10), we obtain the approximation 



exp(-i%(t)) = v^. 



(20) 
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